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Abstract 

We consider an electron-acoustic phonon coupling mechanism associated with the dependence 
of crystal dielectric permittivity on the strain (the so-called Pekar mechanism) in nanostructures 
characterized by strong confining electric fields. The efficiency of Pekar coupling is a function of 
both the absolute value and the spatial distribution of the electric field. It is demonstrated that 
this mechanism exhibits a phonon wavevector dependence similar to that of piezoelectricity and 
must be taken into account for electron transport calculations in an extended field distribution. 
In particular, we analyze the role of Pekar coupling in energy relaxation in silicon inversion layers. 
Comparison with the recent experimental results is provided to illustrate its potential significance. 

PACS numbers: 72.10.-d, 72.10.Di, 73.63.Hs 
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The electron-phonon interaction is one of the fundamental problems in solid-state physics. 
For coupling with acoustic phonons in particular, much attention has been devoted to the 
two main mechanisms in semiconductors: the deformation potential and piezoelectric in- 
teraction. In the presence of an external electric field E, however, an additional process 
appears associated with the dependence of dielectric permittivity on the strain. Due to this 
dependence, an acoustic phonon can induce an effective ac electric field component and sub- 
sequently an additional coupling with electrons. This mechanism was initially introduced by 
p ekar 0, 0] for the problem of sound amplification by drifting electrons. It was shown later 
that the Pekar mechanism of electron-phonon interaction is related directly to the elec- 
trostriction effect. The induced ac electric field has a piezo-like dependence on the phonon 
wavevector q and can roughly be characterized by an effective piezo-constant 

Peff ~ e 2 Ep, (1) 

where p is the photoelasticity constant and e is the dielectric permittivity of the unstrained 
crystal. In ordinary bulk crystals with E < 10 4 V/cm, j3 e ff is quite small. As a result, Pekar 
initially concentrated on the materials with a very large e (e.g., e ~ 2000 for BaTiOs) and 
electrostriction effect. Efficient acousto-electric coupling for such materials was observed 

n 

experimentally in the seventies |^]. 

As the dimension of the sample structure shrinks, interaction mechanisms that are not 
allowed in bulk materials can manifest themselves. For instance, modulation of the quantum 
well width and effective mass by strain gives rise to the so-called macroscopic deformation 
potential |5(. Similarly, the "traditional" mechanisms can exhibit different features. In 
this paper, we demonstrate that the Pekar mechanism is essential for nanostructures with 
strong confining electric fields (f| 3]. The interaction efficiency is dependent on the spatial 
scale of the electric field localization and the phonon wavelength. As a specific example, 
we consider low-temperature electron energy relaxation in n-type (100) Si inversion layers. 
Our calculation shows that due to the Pekar mechanism, the dissipated power in these non- 
piezoelectric materials can exhibit a piezo-like T 3 dependence. This finding provides a clear 
explanation for a recent experimental observation of similar dependence 

Let us start with a general analysis of the Pekar mechanism for the case of a nonuniform 
electric field distribution. The variation of dielectric permittivity under strain Uij can be 
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written in the form 

Seij = efjPiklmUlm, (2) 

where Pikim is the photoelasticity tensor (we could use the electrostriction tensor instead; 
however, the use of photoelasticity tensor allows one to extract explicitly the dependence 
of electrostriction on e, which follows from the Clausius-Mossotti model |H|). The form of 
photoelasticity tensor is determined by the symmetry of the crystal. In the diamond and 
zinc-blende structures, there are only three independent nonzero components of this tensor 
denoted usually as pn, pvz, and P44 j9j. Moreover, = 5ijS. 

The electrostriction effect in a static external field E leads to appearance of an effective 
field E = — V0. The equation for potential <fi can be derived from the equation for the 
electric displacement D: 

V •£> = (), Di = £ij Ejp, (3) 

where E^ = E + E is the total electric field in the crystal. Assuming that the strain is 
small, we obtain 

V^=-^{5e ik E k ). (4) 

In the simplest situation, the external electric field is E — (0, 0, E(z)), which corresponds to 
the confining field in a quantum well. We also assume an elastically uniform medium with the 
functional dependence of Uij,5eij ~ exp[i(—Lat + q z z + qnQ)] and = <f>(z) exp[i(— ut + q\\g)], 
where q\\ = (q x , q y ) and q = (x, y) is the coordinate vector in the quantum well plane. Under 
these conditions, Eq. (J3J) simplifies as 

- -G(z), G(z) = [E(z)^— + 6e zz —^)- } (5) 



dz 2 11 y dxi dz J e 

resulting in 

1 r+OO 1 rz 

<P = -—e q w z J G{z')e- q w z dz' - —e~ q w z J G(z')e^ z dz'. (6) 

It is important to note two important features of 0. First, <fi is not a plane wave as a 
function of z, which is a direct result of the assumed nonuniform character of the electric 
field. Second, strongly depends on the spatial domain of electric field localization, d. For 
q z d ^> 1 and qud 3> 1, the estimate of Eq. ((TJ) applies. For q z d 1 and qud <^ 1, however, 
the induced potential is suppressed substantially with a factor qd. 
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Let us consider a specific example, namely, the process of energy relaxation in n-type 
(100) Si inversion layers at low temperatures. The built-in confining electric fields in such 
structures can be as high as several MV/cm. In this case, the role of Pekar mechanism is 
particularly important since silicon itself is not a piezoelectric material and the Pekar con- 
tribution is expected to be dominant at low temperatures, where the deformation potential 
interaction is less effective. Of course, the Pekar mechanism also provides a contribution to 
the momentum relaxation rate that determines the mobility. However, at low temperatures 
the momentum relaxes mainly due to the elastic scattering by various imperfections, and the 
contribution of phonons is hardly measurable. Thus, we concentrate on energy relaxation. 

The electron potential well near the Si/SiC>2 interface and the electric field associated 
with the inversion layer are shown schematically in Fig. 1. The strong static electric field 
E = —d(j)(z)/dz associated with the inversion layer confines the electrons in a thin silicon 
layer near the interface. The electrostatic potential <f)(z), quantized electron levels and 
wavefunctions ip(r) = y(z) exp(ikp) are determined from the self-consistent solution of 
Poisson and Schrodinger equations p. In the following, we assume that only the ground 
electron subband is populated at low temperatures. 

The induced potential is determined by Eq. (JSJ). The boundary conditions at the Si/SiC>2 
interface (z = 0) under a small strain are: 

= 0(ms) ; £ (s) E (s) + fe(*)£ (s) = e {ins) Ef ls) + 5e^ s) E {ins \ (7) 

where E z = —d<p/dz. For the considered geometry, the following photoelasticity terms are 
relevant: 

5e xz = 2e 2 p AA u xz , 

5e yz = 2e 2 p 44 u yz , (8) 
Se zz = e 2 (p n u zz + Pniuyy + u xx )), 

where x, y, z are the symmetry axis of the crystal. For simplicity, we disregard the mismatch 
of elastic and photoelasticity constants in Si and SiC>2 layers. We also assume that the actual 
phonon wavelength exceeds the thickness of the depletion layer in silicon (see Fig. 1) and yet 
is much less than the thickness of the SiC>2 layer. The former assumption can be justified 
at low enough temperatures, where the contribution of the Pekar mechanism exceeds that 
of the deformation potential. In case when the latter condition is violated, the induced 
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potential becomes weak in accordance with the general analysis given above. In fact, this 
restricts the parameters of actual structures where the Pekar mechanism is important. 
Under these assumptions, / dzx 2 (z)(p(z) « 0(0), leading to 

£ (s) £ (ins) E ( ( q \ \ 

= £ (s) +£ (ins) qil+ \ qz [(P^W + ^aPuK + [pu^ + ipi2 J K<?||) j . (9) 

Here, E s and u s are the electric field and the phonon displacement at the silicon side of the 
interface, respectively. 

Using Eq. JOJ), we calculate the electron scattering rates W^p for electron transition 
k — » k' with absorption (+) or emission (— ) of an acoustic phonon: 

-rrp'ilTP \1 ( F (s) F (ins) \ 2 \ N +1 

W2t = -y^~ [ £ (s )+£{ms) ) I °M I 2 S^8{e k - e kl =F %u q ) ^ (10) 

where is the electron energy and N q represents the Planck population of the phonon modes 
characterized by the lattice temperature T. The index I of the function | ai(8) | denotes the 
type of the phonon modes. In particular, 

I a LA {9) \ 2 = (p 2 n + {Apl 4 - 2p 2 u + 2p nPl2 ) sin 2 9 + {p 2 u + p 2 12 - 2 PnPl2 - Ap 2 M ) sin 4 B) , 

(11) 

\oc TA {0) | 2 = (p44(sm 2 ^-cos 2 ^) 2 + (p 11 -pi 2 ) 2 sin 2 0cos 2 ^), (12) 

where 9 is the angle between the phonon wave vector q and the z axis, V and p are the 
normalizing volume and density, respectively, and uj q = s^ t q is the phonon frequency, olta 
corresponds to the contribution of the transverse phonons with a vertical polarization (u lies 
in the plane formed by the z axis and q); the contribution of the phonons with a horizontal 
polarization (it is parallel to the interface) is zero. The power Q dissipated by electrons is 
given as [ill ]: 

Q = gS- 1 ^^ ~ evWlt + W^kM (1 - f Te (k')) , (13) 

where g is the electron degeneracy (g = 4 in a (100) inversion layer taking into account 
both the spin and valley degeneracy) and S is the cross sectional area for normalization. In 
the expression for Q, we assume that the electron-electron interaction establishes a Fermi 
distribution for electrons with an electron temperature T e . 

At a low temperature, the small-angle (Bloch-Griinizen) scattering regime is realized. For 
degenerate electrons after the standard transformations, we obtain a piezo-like temperature 
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dependence: 



Q = C(T e 3 -T 3 ), 

1 m 2 (E ) 2 e 2 ( r( s )r( ins ) \ 2 
C = ( C LA + C TA) ^^ jm 5^ , (14) 



C L A = \ 77 ( 4 ^44 + 2j?llPl2 + + 3p 
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^he values of the photoelasticity tensor are pn = —0.093, pu = 0.026, 
l^ . As mentioned previously, the electron density and electric field at 



si 8 

Cta = ( 4 P44 + (Pn -Pnf 
s t ° 

where m is the electron effective mass for the in-plane motion, n e is the electron concentration 
in the inversion layer, and ( is the Riemann zeta-function, £(3) ~ 1.20. 

For a numerical estimation, we use the following material parameters: e$i = 12, p = 
2 g/cm 3 , £si0 2 = 4, si = 9 x 10 5 cm/s, s t = 5.4 x 10 5 cm/s, and m = 0.19mo [for the 
(100) interface], 
and p 44 = —0.05 

the interface must be determined by a self-consistent procedure. As a rough estimate, we 
adopt E s = 3 x 10 5 V/cm and n e = 5 x 10 11 cmT 2 : , which provide C = 10" 4 W/K 3 m 2 . 
This corresponds to the energy relaxation time of about 1 ms at the electron temperature 

oi i: :::z:l . Slfan .-^^ ^ t — a, 

the cubic dependence of Q on T e was observed at low temperatures, which was followed by 
a T e 5 dependence for T e > 0.6 K characteristic for the deformation potential coupling. The 
authors of Ref. attributed the T e 3 dependence to the appearance of effective piezoelectric 
properties due to the specific structure of the interface. Indeed, the interface reduces the 
symmetry of the system which can give rise to nonzero piezoelectricity in the vicinity of 
interface. However, for the phonon wavelength exceeding the thickness of this "piezoelectric" 
layer, the piezoelectric potential induced by the phonon is suppressed (this is similar to the 
case of the Pekar potential for narrow regions of electric field localization). Therefore, the 
observed cubic dependence is more likely to be due to the Pekar mechanism. Note also, that 
the obtained above value of C is close to that measured experimentally 0. 

In summary, we show that the Pekar mechanism of electron-phonon interaction (related 
to the electrostriction effect) can become important in nanostructures due to the presence of 
strong confining electric fields and must be considered along with the deformation potential 
and piezoelectric mechanisms. The effectiveness of Pekar coupling depends on both the 



FIG. 1: (a) Schematic illustration of the potential well in a n-type Si inversion layer. The thickness 
of the depletion layer is denoted by dd ep i- (b) z-dependence of the confining electric field in the 
inversion layer. 

absolute value and the spatial distribution of the electric field. An estimation of power 
dissipation by this mechanism in a silicon inversion layer is in good agreement with a recent 
experiment. 
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